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ABSTRACT 

o\ 

. Unparticle physics has been an active field since the seminal work of Georgi. Recently, many 

' constraints on unparticles from various observations have been considered in the literature. In par- 

fSJ , ticular, the cosmological constraints on the unparticle dark component put it in a serious situation. 

In this work, we try to find a way out of this serious situation, by including the possible interaction 
between dark energy and the unparticle dark component. 
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I. INTRODUCTION 

Recently, the so-called unparticle physics has been an active field since the seminal work of Georgi ■ 
It was based on the hypothesis that there could be an exact scale invariant hidden sector resisted at a 
high energy scale. A prototype model of such a sector is given by the Banks-Zaks theory which flows 
to an infrared fixed point at a lower energy scale Au through dimensional transmutation Recently, 
there has been a flood of papers on the unparticle phenomenology. We refer to e.g. [1, 01 for some brief 
reviews. In fact, soon after the seminal work of Georgi [l], l^, some authors considered the constraints 
on unparticles from various observations, such as the collider experiments the new long range force 
experiments 01 , the solar and reactor neutrinos data [1] , as well as the observations from astrophysics 
and cosmology 13, H, Hi [0, 113 ■ These observations put fairly stringent constraints on the unparticle. 

The unparticle does not have a definite mass and instead has a continuous spectral density as a 
consequence of scale invariance [l|,[3] (see also e.g. [Ulillll), 

p{P^) = AaJ{P°)e{P^){P^Y--\ (1) 

where P is the 4- momentum; is the normalization factor; du is the scaling dimension. When du = 1, 
Eq. reduces to the familiar one of a massless particle. The scaling dimension c?u is constrained 
by the unitarity of the conformal algebra and the requirement to avoid the singular behavior [l|, Q]. 
The theoretical bounds are 1 < du < 2 (for bosonic unparticles) or 3/2 < du < 5/2 (for fermionic 
unparticles) Q. Typically, values 1 < du < 2 have been extensively considered in the literature. 
The pressure and energy density of thermal unparticles have been derived in They are 

T c(du) 



A, = (2d„ + i)<,.r« eg^, (3) 

where C (du) = B (3/2, du) T (2du + 2) C (2du + 2), while B, T, C are the Beta, Gamma and Zeta func- 
tions, respectively. These are the results for the bosonic unparticles. The pressure and energy density 
of fermionic unparticles can be obtained by replacing C (du) by [l - 2-(2'iu+i)] q (du). Therefore, the 
equation-of-state parameter (EoS) of both the bosonic and fermionic unparticles is given by 



Pu 1 



Pu 2du + 1 ' 



(4) 



When du = 1, we have Wu = 1/3, which is the same as radiation. When du oo, we find Wu 0, which 
approaches the EoS of pressureless matter. In the intermediate case, the EoS of unparticles is different 
from the one of radiation or cold dark matter, and generically lies in between. Since the unparticle 
interacts weakly with standard model particles, it is "dark" in this sense. Some authors regard the 
unparticle as a candidate of dark matter [l^ (see also [111, [l^l)- Noting that the EoS of unparticles 
(wu > 0) is different from the one of cold dark matter (wm =0), we instead call it "dark component" to 
avoid confusion. 

In the cosmological constraints on the unparticle dark component have been considered, by using 
the type la supernovae (SNIa), the shift parameter of the cosmic microwave background (CMB), and the 
baryon acoustic oscillation (BAO). The authors of [ll| found that du > 60 at 95% confidence level (C.L.) 
for the AUDM model in which the unparticle is the sole dark matter. As mentioned above, however, 
the theoretical bounds on unparticles are 1 < du < 2 (for bosonic unparticles) or 3/2 < du < 5/2 (for 
fermionic unparticles) [3]. The situation is serious. Even for the AUCDM model in which the unparticle 
dark component co-exists with cold dark matter, they found that the unparticle dark component can at 
most make up a few percent of the total cosmic energy density if du < 10; so that it cannot be a major 
component In fact, it is easy to understand these results. Since the unparticle dark component scales 
as a-^Cn-iiiu) whereas cold dark matter scales as (here a is the scale factor of the universe), for > 0, 
the energy density of the unparticle dark component decreases faster than the one of cold dark matter. 
So, it is not surprising that the energy density of the unparticle dark component cannot be comparable 
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with cold dark matter at the present epoch. To make a considerable contribution to the total cosmic 
energy density, should be very large to make ~ so that the unparticle dark component could 
mimic cold dark matter. 

In this work, we try to find a way out of the serious situation mentioned above. Our physical motivation 
is very simple. If dark energy, the major component of the universe, can decay into unparticles, the energy 
density of the unparticle dark component should decrease more slowly. Since both dark energy and 
unparticle are unseen, the interaction between them is not prevented. If the dilution rates of unparticle 
dark component and cold dark matter are comparable, it is possible to have a considerable energy density 
of the unparticle dark component at the present epoch, without requiring to be very large. In fact, 
this is similar to the key point of the interacting dark energy models which are extensively considered in 
the literature to alleviate the cosmological coincidence problem (see e.g. [3 HI] and references therein). 

Here, we consider a flat universe which contains dark energy, the unparticle dark component and 
pressureless matter (including cold dark matter and baryons). We assume that dark energy and unparticle 
dark component exchange energy according to (see e.g. [25| and references therein) 

Pdc + 3iJpde(l + Wdc) = -Q, (5) 
Pu + 3i//9u(l+Wu) =Q, (6) 

whereas pressureless matter (including cold dark matter and baryons) evolves independently, i.e., 

p„, + 3Hp„, = 0. (7) 

So, the total energy conservation equation ptot + 3i7/9tot(l + Wtot) = is preserved. H = a/a is the 
Hubble parameter; a = (1 + z)^^ is the scale factor of the universe (we set oq = 1); the subscript "0" 
indicates the present value of the corresponding quantity; z is the redshift; a dot denotes the derivative 
with respect to cosmic time t. The interaction forms extensively considered in the literature (see for 
instance [l3.[25j| and references therein) are Q oc Hp^ji, Hp^^, -ffptot, i^Pm<j> (where = SttG), and so on. 
In this work, for simplicity, we consider the interaction term 

Q = 3aHp^, (8) 

where a is a constant. So, Eq. ® becomes pu + iHp^ (l + w^^) = 0, where w'^ = w-^ — a. It is easy to 
find that 



-3(i+<0. (9) 



On the other hand, from Eq. ([7]), we have pm = PmO<i~^- For convenience, we introduce the fractional 
energy density Vli = {SirGpi) / (3H^), where i = de, u and m. 

In the following sections, similar to [ll| . we consider the cosmological constraints on the unparticle 
dark component which interacts with dark energy, by using the observations of SNIa, the shift parameter 
R of CMB, and the distance parameter A from BAO. In Scc.|TTl we present the cosmological data and the 
methodology used to constrain the models. In Sec. IIIII and llVi we consider the observational constraints 
on the interacting AUCDM model and the interacting XUCDM model, respectively. As expected, we find 
that the cosmological constraints on the unparticle dark component can be significantly relaxed, thanks 
to the possible interaction between dark energy and unparticle dark component. The tension between 
theoretical bounds and cosmological constraints could be removed. Finally, some concluding remarks are 
given in Sec. |Vl 

II. COSMOLOGICAL DATA 

In this work, we perform a analysis to obtain the constraints on the model parameters. The data 
points of the 307 Union SNIa compiled in [l^ are given in terms of the distance modulus Pobs{zi)- On 
the other hand, the theoretical distance modulus is defined as 



Pth{zi) = Slogio DL{zi) + po, 



(10) 
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where hq = 42.38 — Slogj^g h and h is the Hubble constant Hq in units of 100 km/s/Mpc, whereas 
in which E = H/Hq; p denotes the model parameters. The from the 307 Union SNIa are given by 

where a is the corresponding Ict error. The parameter /io is a nuisance parameter but it is independent of 
the data points. One can perform an uniform marginalization over fiQ. However, there is an alternative 
way. Following [l6l . [l7j , the minimization with respect to fj,o can be made by expanding the xfj, oi Eq. (|12p 
with respect to fiQ as 

xl{p) = A-2fioB + filC, (13) 

where 

[^J'obs{zi) - fithizi- fj.0 = o,p)f 



Mp) = E ■ 



^ f^obsizt) - lJ^th{zi;fio = 0,p) \r 1 

""^^^^^ <J^) ' ^ = ^<J^)- 

Eq. (|13p has a minimum for /io = -B/C at 

X^(p) = i(p)-^. (14) 

Since xfi min — xfi mm obviously, we can instead minimize which is independent of /ip. The shift 
parameter R of CMB is defined by [H 

H = nl^f^. (15, 

where the redshift of recombination = 1090 which has been updated in the Wilkinson Microwave 
Anisotropy Probe 5-year (WMAP5) data ^20j. The shift parameter R relates the angular diameter 
distance to the last scattering surface, the comoving size of the sound horizon at and the angular 
scale of the first acoustic peak in CMB power spectrum of temperature fluctuations [3, [HI . The value 
of R has been updated to 1.710 ± 0.019 from the WMAP5 data [2^. On the other hand, the distance 
parameter A of the measurement of the BAO peak in the distribution of SDSS luminous red galaxies [2l| 
is given by 

where Zb = 0.35. In [2^1, the value of A has been determined to be 0.469 (ns/0.98)"°-35 ± 0.017. Here the 
scalar spectral index is taken to be 0.960, which has been updated from the WMAP5 data [13] ■ So, 
the total is given by 

= XI + XCMB + X%AO^ (17) 

where xf, is given in Eq. Xcmb = {R- RobsY /<Jr and x%ao = (A- Aobsf /<t\- The best-fit model 
parameters are determined by minimizing the total x^- As is well known, the 68% C.L. is determined by 
Ax^ = X^ - xlitn < 1-00, 2.30, 3.53, 4.72 and 5.89 for n = 1, 2, 3, 4 and 5, respectively, where n is the 
number of free model parameters. Similarly, the 95% C.L. is determined by Ax^ = — Xmin — 4.00, 
6.17, 8.02, 9.70 and 11.3 for n = 1, 2, 3, 4 and 5, respectively. 

In the following sections, we consider the observational constraints on the interacting AUCDM model 
and the interacting XUCDM model, respectively, by using the methodology presented here. 



1 

_Zb Jo 



dz 

W). 



(16) 
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III. COSMOLOGICAL CONSTRAINTS ON THE INTERACTING AUCDM MODEL 

In this section, we consider the interacting AUCDM model in which the role of dark energy is played by 
the decaying cosmological constant (vacuum energy), as in the well-known A(i)CDM model ^26,. .27.. i28j. 
In this case, Eq. ([SJ becomes 

PA = -iaHp^ = -3ap^oHa-^(^+<"). (18) 

Noting that the case without interaction (a = 0) was considered in llj, there are two other different 
cases with w^^ 7^ —1 and w'^^ = — 1 for this differential equation. We will study them one by one. 




FIG. 1: The 68% C.L. (black solid lines) and 95% C.L. (red dashed lines) contours in the fluo — du plane and the 
Q — du plane for the interacting AUCDM model with w'^ 7^—1- The best fit is indicated by the red solid point. 



A. The interacting AUCDM model with w'f / -1 



If wf ^ 



-1, the solution for Eq. (|18p is given by 

«PuO ^_3(l+<«)^^ 



PA = 



1 



„off 



const., 



(19) 



where const, is an integration constant. Obviously, this solution diverges for — —1 which will be con- 
sidered in the next subsection. Inserting Eq. into the Friedmann equation = {pA + Pu + Pm) /3 
and requiring E(z = 0) = 1, one can determine the integration constant. Finally, wc find that 



f^mo(l + z)'+ 1 + 



1 



„cff 



1-a, 



1 



„cff 



(20) 



There are 4 free model parameters, namely, rimOj ^uo, du and a. By minimizing the corresponding total 
in Eq. we find the best-fit parameters f},„o = 0.280, Q^o = 0.005, = 1.722 and a = 0.116, 

while Xmin — 311.924. In Fig.[Tl we present the corresponding 68% and 95% C.L. contours in the iluo — du 
plane and the a — du plane for the interacting AUCDM model with 7^ — 1 , while the other parameters 
are taken to be the corresponding best-fit values. Obviously, from the left panel of Fig. [H we see that 
the unparticle dark component can at most make up a few percent of the total cosmic energy density, so 
that it cannot be a major component in this case. 
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B. The interacting AUCDM model with wf = -1 

As mentioned above, the solution Eq. ITOl) diverges for w'^ = — 1. Therefore, one is required to consider 
this case separately. If w"^ = —1, we find that 

+ + (21) 

So, a is no longer a free model parameter. In this case, Eq. ^ becomes pu = PuO- Now, the solution for 
Eq. is given by 

PA = — SapuO In a + const., (22) 



where const, is an integration constant. Then, inserting Eq. ()22p into the Friedmann equation and 
requiring E{z = 0) = 1, one can determine the integration constant. Finally, we find that 

= r!,,o(l + z f + Sriuo (^1 + + + (1 - r!mo) • (23) 

In this case, there are 3 free model parameters, namely, rimOi ^uo and d^. By minimizing the correspond- 
ing total in Eq- ([T7)) . we find the best-fit parameters ilmo = 0.270, fluo = 0.020 and = 0.273, while 
Xmin — 313.168. In Fig. ^ we present the corresponding 68% and 95% C.L. contours in the iluo — du 
plane for the interacting AUCDM model with w'^^ = —1, while the other parameters are taken to be 
the corresponding best-fit values. It is a pleasure to find that the unparticle dark component can have 
a considerable contribution to the total cosmic energy density {fluo > 0.1) in the 95% C.L. region while 
1 < c^u < 2_('for bosonic unparticles) or 3/2 < du < 5/2 (for fermionic unparticlcs). The serious situation 
found in \1m could be relaxed in this case. 



14 
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FIG. 2: The 68% C.L. (black solid line) and 95% C.L. (red dashed line) contours in the fluo — du plane for the 
interacting AUCDM model with w"^ = —1. The best fit is indicated by the red solid point. 
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IV. COSMOLOGICAL CONSTRAINTS ON THE (INTERACTING) XUCDM MODEL 



In this section, we consider the (interacting) XUCDM model in which the EoS of dark energy Wx is a 
constant. We can recast Eq. (O as 

t^ + -p.(l + u..) = -%.oa-(^+<"). (24) 
da a a 

In the following, we first consider the case without interaction (a = 0). Then, we consider two different 
cases with w'^^ ^ and w"^ = Wx for the differential equation 



A. The XUCDM model without interaction 



It is natural to see first what will happen if there is no interaction between dark energy and unparticle 
dark component. In this case, a — 0. One can easily write down 

E'^ = n^i + zf + n^oii + z)3(i+-") + (1 - nn,o - ^uo) (i + zf'^^+^^l (25) 

There are 4 free model parameters, namely, flmo, f^uOi and w^- By minimizing the corresponding total 
in Eq. p?)) . we find the best-fit parameters Qnio = 0.283, iluo ~ 0.0, ~ 0.540 and = —1.006, 
while Xmin — 311.981. In Fig. [31 we present the corresponding 68% and 95% C.L. contours in the QuQ — du 
plane and the — du plane for the XUCDM model without interaction, while the other parameters are 
taken to be the corresponding best-fit values. Obviously, from the left panel of Fig. [H we see that the 
unparticle dark component can at most make up a few percent of the total cosmic energy density, so that 
it cannot be a major component in the case without interaction. 




FIG. 3: The 68% C.L. (black solid lines) and 95% C.L. (red dashed lines) contours in the fluo — du plane and the 
Wx — du plane for the XUCDM model without interaction. The best fit is indicated by the red solid point. 



B. The interacting XUCDM model with wf / Wx 



If w° ^ Wx, the solution for Eq. ([24|l is given by 



^ ^^"0 a-3(i+<") + const, x a-3(i+-^) , (26) 

10™ - Wyr 
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where const, is an integration constant. Obviously, this solution diverges for w'^ = which will be con- 
sidered in the next subsection. Inserting Eq. into the Friedmann equation = {p^ + Pu + Pm) /3 
and requiring E(z = 0) = 1, one can determine the integration constant. Finally, we find that 



a 



1 - i^n,0 - II 



(27) 



There are 5 free model parameters, namely, flmo, f^uO: f^u, « and w^- By minimizing the corresponding 
total in Eq. HT]), we find the best-fit parameters Q^o = 0.277, QuO = 0.214, = 13.478, a = 0.589, 



and 



— —1.666, while 



xl^^^ = 311.077. In Fig. HI we present the corresponding 68% and 95% C.L. 
contours in the fluo — plane and the a — du plane for the interacting XUCDM model with w'^ ^ w^, 
while the other parameters are taken to be the corresponding best-fit values. It is fortunate to find 
that the unparticle dark component can have a sizable contribution to the total cosmic energy density 
(0.1 < n^o < 0.2 in the 68% C.L. region and 0.1 < Qno < 0.25 in the 95% C.L. region) while 1 < du < 2 
(for bosonic unparticles) or 3/2 < d^ < 5/2 (for fermionic unparticles) . The serious situation found 
in [ll[ could be significantly relaxed in this case. 
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FIG. 4: The 68% C.L. (black solid lines) and 95% C.L. (red dashed lines) contours in the Cluo — du plane and the 
Q — du plane for the interacting XUCDM model with w'^ 7^ lUx. The best fit is indicated by the red solid point. 



C. The interacting XUCDM model with wf = 

As mentioned above, the solution Eq. (|26p diverges for w'^ = w^- Therefore, one is required to consider 
this case separately. If w'^ = Wx, we find that 

a = Wu-Wx= ^ - Wx- (28) 

So, a is no longer a free model parameter. In this case, the solution for Eq. is given by 

Px = {const. - 3apno In a) a-3(i+«'x) ^ (29) 
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where const, is an integration constant. Inserting Eq. (|29p into the Friedniann equation and requiring 
E(z = 0) = 1, one can determine the integration constant. FinaUy, we find that 

There are 4 free model parameters, namely, ftmo, ^uo, and w^. By minimizing the corresponding total 
in Eq. dTT]), we find the best-fit parameters f^mo = 0.270, QuO = 0.0, du = 1.969 and = -0.954, 
while Xmin — 313.138. In Fig. we present the corresponding 68% and 95% C.L. contours in the 
^^uo ~ dn plane and the Wx — rfu plane for the interacting XUCDM model with = Wx, while the other 
parameters are taken to be the corresponding best-fit values. It is a pleasure to find that the unparticle 
dark component can have a considerable contribution to the total cosmic energy density (fiuo > 0.1) 
in the 95% C.L. region while 1 < du < 2 (for bosonic unparticles) or 3/2 < d^ < 5/2 (for fermionic 
unparticles). The serious situation found in [ll| could also be relaxed in this case. 
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FIG. 5: The 68% C.L. (black solid lines) and 95% C.L. (red dashed lines) contours in the fluo — du plane and the 
Wx — du plane for the interacting XUCDM model with The best fit is indicated by the red solid point. 



(1 — ^mo) + 3f7uO 



2d„ 



Wx ln(l + z) 



(1 + ^) 



3(l+u 



(30) 




V. CONCLUDING REMARKS 

Unparticle physics has been an active field since the seminal work of Georgi [H, g. Recently, many 
constraints on unparticles from various observations were considered in the literature. In particular, the 
cosmological constraints on the unparticle dark component put it in a serious situation [ll[ . In the present 
work, we try to find a way out of this serious situation, by including the possible interaction between 
dark energy and the unparticle dark component. In fact, as we have shown above, the cosmological 
constraints on the unparticle dark component can be relaxed in the interacting AUCDM model with 
w'^ = — 1, and both the interacting XUCDM models with w'^ ^ Wx and w'^ = Wx- In these interacting 
models, the unparticle dark component can have a considerable contribution to the total cosmic energy 
density (namely fl^o > 0.1; even fl^o can be larger than 0.2 in the interacting XUCDM model with 
w'^ 7^ Wx) while 1 < du < 2 (for bosonic unparticles) or 3/2 < du < 5/2 (for fermionic unparticles), 
thanks to the possible interaction between dark energy and unparticle dark component. The tension 
between theoretical bounds and cosmological constraints could be removed. 
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It is easy to understand these results physically. When dark energy, the major component of the 
universe, can decay into unparticles, the energy density of the unparticle dark component should decrease 
more slowly. Since both dark energy and unparticle are unseen, the interaction between them is not 
prevented. If the dilution rates of unparticle dark component and cold dark matter are comparable, it 
is possible to have a considerable energy density of the unparticle dark component at the present epoch, 
without requiring to be very large. In fact, this is similar to the key point of the interacting dark 
energy models which are extensively considered in the literature to alleviate the cosmological coincidence 
problem (see e.g. [13, Hll and references therein). 

As shown above, the cosmological constraints on the unparticle dark component can be relaxed in the 
interacting AUCDM model with w'^^ = — 1, and both the interacting XUCDM models with w'^ ^ 
and = Wx- However, one might note that the unparticle dark component with cosmological constant 
is constrained by the observational data relatively more severe than in the case with dark energy. This 
is also easy to understand. Notice that the EoS of dark energy is a free model parameter, whereas 
the EoS of the cosmological constant is a constant —1. Therefore, the number of free model parameters 
in the interacting XUCDM model is larger than the one in the interacting AUCDM model. As is well 
known, statistically, the constraints should be looser in this case. 
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